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Abstract 

A review of earlier fluid dynamical calculations with QGP show a softening of the 
directed flow while with hadronic matter this effect is absent. The effect shows up 
in the reaction plane as enhanced emission which is orthogonal to the directed flow. 
Thus, it is not shadowed by the deflected projectile and target. As both of these 
flow components are in the reaction plane these form an enhanced 'elliptic flow' 
pattern. Recent experimental data at 11 AGeV and above show the same softening, 
hinting at QGP formation. 



1 Introduction 



It is widely known that fluid dynamics (FD) is governed by the Equation 
of State (EoS) of the matter, and the analysis of the resulting flow patterns 
turned out to be one of the best tools to extract the EoS from the outcome of 
a heavy ion collision. The final event shape must carry the information about 
the pressure development during the collision including the early stages of the 
collision. 



The phase transition to the Quark-Gluon Plasma (QGP) is connected to a de- 
crease in pressure according to most theoretical estimates, not only in strong 
first order phase transition models, but even if we have a smooth but rapid 
gradual transition. This reduced pressure and temperature around the phase 
transition threshold is known for long (see e.g. [1]), and it is emphasized re- 
cently as a QGP signal, i.e., the "soft point" of the EoS, which should be 
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possible to observe in excitation functions of collective flow data[5,6]. It was 
pointed out earlier [3] that the decrease in the out of plane (squeeze out) emis- 
sion is even more sensitive to the pressure drop and it decreases due to plasma 
formation stronger than the inplane collective flow. 

Here we want to point out another consequence of the same softening in the 
EoS, which is a new distinct flow pattern, which can be seen in almost all 
theoretical FD calculations with QGP formation, but which was not discussed 
earher[2-5,7]. 



2 Third flow component in fluid dynamics 

Up to now two basic flow patterns are predicted and detected: (i) the directed 
transverse flow in the reaction plane, or side-splash, or bounce off, which is 
most frequently presented on the well known P^. vs. y diagram and seen at 
all energies in heavy ion collisions from energies of 30 A.MeV to 165 A.GeV, 
[10,13,11,12] and the (ii) squeeze-out effect which is an enhanced emission of 
particles transverse to the reaction plane at center of mass (CM) rapidities. 

At lower energies the directed transverse flow resulted in a smooth, linear 
vs. y dependence at CM rapidities. This straight line behavior connecting the 
maximum at yproj and the minimum at ytarg was so typical that it was used 
to compare flow data at different beam energies and masses. 

If QGP is formed, strong and rapid equilibration and stopping takes place, and 
close to one-fluid behavior is established. Stopping is stronger than expected, 
and Landau's fluid dynamical model is becoming applicable for central colli- 
sions of massive heavy ions. The soft and compressible QGP forms a rather 
flat disk orthogonal to the beam axis which is at rest in the CM system. Then 
this disk starts to expand rapidly in the direction of the largest pressure gra- 
dient, i.e., forward and backward. Thus, the not fully Gaussian shape of the 
measured rapidity spectra can be interpreted as a fluid dynamical bounce back 
effect (Landau model) in contrast to the transparency otherwise assumed in 
kinetic models. Unfortunately in central collisions we can not distinguish the 
two effects from one - another. Both lead to a spectrum elongated in the beam 
direction. 

At small but finite impact parameters, however, this disk is tilted, and the 
direction of fastest expansion, will deviate from the beam axis, will stay in the 
reaction plane, but point in directions opposite to the standard directed trans- 
verse flow. Since pressure does not play a role in transparency, transparency 
cannot explain such deviation from the beam direction! This third flow com- 
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ponent develops purely from the large pressTire gradient at ftill stopping of 
the strongly Lorentz contracted intermediate state. So, at the same time as 
the primary directed flow is weakened by the stronger Lorentz contraction 
at higher energies, this third flow component is strengthened by increased 
Lorentz contraction. These two flow components together form the 'elliptic 
flow'. [11, 16,17] 

On the Px vs. y diagram [14] this component shows up as a smaller, negative 
flow component at small CM rapidities. Such a third flow component is seen 
clearly in Fig. 3 of [3] (see Fig. 1, lower part), Fig. 8 of [4], Fig. 6b of [5] 
and Fig. 6 of [7] at or slightly below 0.5 y/Ucm if QGP formation was allowed 
during the calculation. In sharp contrast, the solutions with hadronic EoS 
did not show this effect, and the maximum and minimum of the Px curve 
could be connected with a rather straight hne. This straight line behavior is 
typical at all flow results below 11 A.GeV beam energy (Fig. 2). In some of 
the FD calculations with QGP the secondary peak at small CM rapidities 
is not seen, but a tendency is obvious, and the deviation from the hadronic 
smooth line behavior is apparent. This can be seen clearly in Fig. 3 of [2], and 
Figs. 6a and 6c of [5]. This indicates that the strength of this effect is also 
impact parameter and beam energy dependent, and the third flow component 
shows a relative maximum at the same energy when the primary directed flow 
is at its minimum [5]. Note that all these FD calculations were done much 
before the experiments. The first quantitative flow predictions [2] preceded 
the experiments by as much as 6 years (!) and gave rather good agreements 
with the data. 

To have a quantitative measure of the softening at small CM rapidities {ycM — 
0) for a symmetric (A+A) collision, on the \Px-, y\ plane we draw the steepest 
straight line, ay, through the CM point which is tangent to the Pxi^U) curve 
(Fig. 1 upper part). The Px{y) curve and the straight line are usually tangent 
to each other at y = 0, and at finite y the Px{y) curve is under the straight 
line, I ay I > \Px{ii)\. In case of softening at low rapidities the straight line is 
tangent to the Pxiy) curve at two points, yi and —yi and at smaller rapidities, 
and at smaller but finite rapidities, < \y\ < \yi\, the Px{y) curve is under 
the straight line, \ay\ > \Px{y)\. The relative deviation of these two quantities, 
S{y) = \ay — Px{y)\/\ay\ has a maximum at some rapidity, < \y = ymax| < 
\yi\, and we will use this quantity, S{y„^a_^), to characterize the softening of 
EoS. If 5" = 100% this means that Px{y) vanishes somewhere between and 
|yi|. If 5" > 100% this means that Px is inverted (shows negative flow) at low 
rapidities. 



Thus in FD calculations, this effect is clearly predicted since 1991, in all cal- 
culations, without being noticed or discussed up to now. Nevertheless, all FD 
calculations are consistent in predicting an observable softening if and only if 
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Publication 



QGP HM 



[5] Au+Au 3.5 A.GeV b = 3fm 

[5] Au+Au 5 A.GeV b = 3fm (soft pt.) 

[5] Au+Au 11.7 A.GeV b = 3fm 

[3] Au+Au 11.6 A.GeV b = 5fm 

[4] Au+Au 11.6 A.GeV 6 = 0- OM^^ 

[7] Au+Au 11 A.GeV 6 = 3fm 6 fm/c 

[7] Au+Au 11 A.GeV b = 3fm 7.2 fm/c 

[2] Pb+Pb 160 A.GeV b = 4fm 3.8 fm/c 



180 ±25% 30 ±30% 



100 ±10% 



160±30% 



50 ±35% 



80 ±10% 50 ±45% 



70 ±30% 



70 ±10% 30 ±30% 



40% 



±20% 



0% 



0% 



Table 1 



Softening in fluid dynamical model calculations which included the possibility of 
QGP formation in the EoS assumed. When Hadronic Matter (HM) EoS was also 
presented in the same calculation the resulting softening is also shown for this case. 
All HM calculations are consistent with zero softening, while calculations with QGP 
show softening values from 50 to 180%. (In [7] the non-equilibrium 3- fluid hydro is 
applied with complex EoS, i.e. hadron gas EoS for the two baryon fluids (no QGP), 
and EoS with QGP-hadron phase transition, which is used only for the third, baryon- 
free fluid. Although the effects like directed flow are provided by baryons the flow 
effect arises from the pressure of the central baryon free component. Nevertheless, 
non-equilibrium effects in 3-fluid hydrodynamical models could mimic the effect 
caused by the QGP in ideal 1-fluid models, and so the comparison of this model 
with other fluid dynamical models is somewhat arbitrary.) 

QGP EoS is included in the calculation. 



3 Deviation from straight line behavior in string models 

Some of the recent flow data are analyzed by the Fourier expansion method 
[15,16]. Unfortunately, it is not possible to convert or compare these data to 
the earlier P^. analysis, but the results indicate that the Vi vs. y/Ucm has a 
straight line behavior at low rapidities in the absence of QGP (or softening) 
in the EoS. Thus, just as in the Px analysis we consider that deviation from 
this straight line behavior indicates a stronger contraction and softening of 
the EoS. To make the comparison more quantitative would only be possible 
if data showing such deviation from the straight line behavior would be an- 
alyzed by both methods simultaneously. Nevertheless, heavy ion collisions at 
lower energies, around 2 A.GeV show straight line behavior according to both 
analyses[19]. Furthermore, small relative deviations from the straight hue be- 
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havior (or softening) are expected to be identical in the two analyses [10] since 



We do not intend to compare FD models and non-equilibrium models like 
string models as theoretically they are based on different fundations. However, 
wc review how these models describe the data. In the last decade string models 
showed a surprising level of resilience by being able to reproduce more or less 
all measured data. Unlike with FD models this always happened after the 
data became known. Of course during this time the models were changed 
essentially and many new ingredients were built in. Some of these changes are 
in fact rather essential, and change the physical picture described by these 
models drastically. One of the most important changes is the introduction of 
string ropes (fused strings, or quark clusters), and only this made it possible to 
describe the formation of massive formed secondaries like strange antibaryons. 
Recent flow data are also stretching the string models flexibility to their hmits, 
while FD models as usual overestimate collective flow data. 

In Table 2 we present a collection of string model results regarding the low 
rapidity softening. 



In an attempt to fit all data with string models in ref. [8] RQMD was particu- 
larly tested to which extent it is able to reproduce NA49 flow data. Fig. 3 [8] 
presents the Vi vs. ij plot for 158 A.GeV Pb+Pb data. The proton data show 
a very strong deviation from the straight line behavior (indicating strong EoS 
softening based on analogies with results). The best flts with RQMD v 2.3 
also reproduce some of the deviation from the straight line behavior, but at 
low rapidities the model result lags behind the SPS data appreciably, while 
it yields also some softening at 2 A.GeV where experiments show no sign of 
such behavior. 

Venus and to some extent RQMD might include some of the genuine flow 
effects as the increasing Lorentz contraction, increasing pressure, and increased 
stopping, but these do not quite reach the extent one can have in an FD model 
with QGP. 

Although, the authors of ref. [8] conclude that in RQMD the deviation from 
the straight line behavior is caused mainly by shadowing, it is difficult to 
observe any systematic viewing all published RQMD results. The shadowing 
effect is also confirmed by an earlier work [9] which actually went more into 
details, and showed that for heavy systems hke Pb+Pb the shadowing starts to 
cause deviations from the straight line behavior only in peripheral reactions 
(6 > 0.756max)- For lighter systems the deviation may show up already at 
smaller impact parameters also. This coincides with the result shown in Fig. 
2 of [8], for Pb+Pb and S+S at 158 A.GeV. 
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Publication 



S 



[18] Venus 4.12 Pb+Pb 158 A.GeV vi 

[12] Venus 4.12 Pb+Pb 158 A.GeV b = 5- 12fm 

[8] RQMD 2.3 Pb+Pb 2 A.GeV 6 = 5 - 8fm Fig. 2 vi 

[10] RQMD 2.3 Au+Au 11.6 A.GeV 

[10] RQMD 1.08 Au+Au 11.6 A.GeV P^ 

[12] RQMD 2.3 Pb+Pb 158 A.GeV b = 8- lOfm P^ 

[8] RQMD 2.3 Pb+Pb 158 A.GeV Fig. 1 vi 

[8] RQMD 2.3 Pb+Pb 158 A.GeV 6 = 5 - 8fm Fig. 2 vi 

[2] QGSM Pb+Pb 158 A.GeV b = 4fm P^ 

[3] QGSM Au+Au 11.6 A.GeV b = 3fm P^ 

[9] QGSM Au+Au 11.6 A.GcV 6 = 1, lOfm P^ 

[9] QGSM Au+Au 11.6 A.GeV b = llfm P^ 

[9] QGSM Au+Au 11.6 A.GeV b = 12fm P^ 



110+15% 



134+15% 



32+20% 



0+25% 



0+20% 



0+30% 



0+10% 



40% 



35% 



50% 



72% 



57% 



0% 



Table 2 



Softening in string model calculations. All QGSM calculations are consistent with 
zero softening, with the exception of extreme peripheral collisions, where the effects 
of shadowing are observable. Both Venus 4.12 calculations yield consistently large 
softening at SPS energies. RQMD results show varying results from zero to 72%. 

From these string model results we can conclude that shadowing effects in 
reality may also result in deviations from the straight line behavior, but in 
heavy colliding systems this happens only in rather peripheral reactions, which 
can be excluded experimentally. In the Pb+Pb case the shadowing effect below 
the used impact parameter cut of 6 < 8fm is rather weak (cf . [9] ) thus cannot 
be solely responsible for the deviation from the straight line behavior. We have 
to add that the previously discussed purely fluid dynamical effects also show up 
in string model results in a weaker form. As this is caused by the considerable 
softening of the EoS due to string ropes and other massive objects, we can 
actually observe some level of convergence between FD- and string models. 
Nevertheless, earlier the authors of string model calculations did not observe 
the possible fluid dynamical component in the explanation. In a recent version 
of RQMD the measured flow pattern could be reproduced by including fluid 
dynamical mechanisms and an EoS. [17] 

In view of the extreme variety of string model results, wc cannot reach any 
conclusion whether string models consequently yield low energy softening or 
not. Furthermore, it is not clear what is the general, physical reason of the 
softening in those string model calculations where such a softening is present. 
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4 Experimental status 



In connection with the string models we discussed already NA49 data. Other 
experiments starting at about 2-11 A.GeV (Figs. 2 and 3) show signs of de- 
viation from the straight line behavior. At 11 A.GeV and particularly at SPS 
energies the vs. y (Fig. 2 bottom) and the v\ vs. y (Fig. 3 bottom) becomes 
almost flat at or close to CM rapidities. Although the inversion of proton flow 
is not observed as predicted by some of the FD calculations, the tendency 
from the straight line behavior is consistently seen in all experiments (Table 
3). 



Publication 


S 


[13] Au+Au 4 A.GeV 


40+50% 


[13] Au+Au 6 A.GcV P^ 


10+35% 


[10] Au+Au 11 A.GeV P^ protons 


30+20% 


[12] Pb+Pb 158 A.GeV P^ semicentr. 


60+50% 


[11] Pb+Pb 158 A.GeV vi 


110+60% 



Table 3 

Softening in directed transverse flow data. With increasing beam energy the ten- 
dency of an increasing softening is observable. At SPS energies the softening becomes 
significant. This general trend coincides with the fluid dynamical predictions with 
QGP EoS. 

The quantitative study of this effect in finer impact parameter resolution may 
even pin down more pronounced deviations in the central rapidity region which 
is important if we want to use flow data to reconstruct the properties of the 
underlying EoS. 



5 Conclusion 

The observed flow patterns at ultra-relativistic energies indicate a new flow 
pattern, which may arise from a highly Lorentz contracted and compressed 
intermediate state as a consequence of an extremely soft EoS. The presently 



7 



observed deviations from the straight hnc behavior go beyond hadronic or 
string model predictions, indicating a softer and more compressed initial or 
intermediate state than hadronic models can accommodate. 

A review of earlier FD calculations which included QGP shows that such an 
effect was present already in all of them to a smaller or larger extent. This 
sign of excessive softening of the EoS may indicate that a larger portion of the 
matter is transformed into a soft phase than assumed in string models. The 
effect shows up in the reaction plane as enhanced emission which is orthogonal 
to the directed flow. Thus, it is not shadowed by the deflected projectile and 
target. [3,9] As both of these flow components are in the reaction plane these 
form an enhanced 'elliptic flow' pattern. 

A review of string model results does not provide us with any clear conclusion 
about the causes or extent of the softening, although the effect is present in 
some of the string model calculations. 

In order to be more quantitative further and more detailed experiments are 
needed, as well as further FD calculations would be beneficial. Fluctuations, 
initial pre-equilibrium transparency, viscosity, etc. may decrease the effect in 
FD calculations also, thus approaching the real reaction mechanism better 
than recent calculations. 
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-1 1 



y/y^ 



Fig. 1. Upper part: Definition of the measure softening, S, describing the devia- 
tion of Px{y) or vi{y) from the straight hne behavior, ay, around CM. S is defined 
as \ay — Px{y)\/\ay\. The lower figure shows a typical example for fluid dynami- 
cal calculations with Hadronic and QGP EoS [3]. QGP leads to strong softening, 
~ 100%. 
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Fig. 2. Results of transverse flow analyses, Px vs. y/ycm, [8,8,10,12]. 158 AGeV 
results show strong softening. 
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Fig. 3. Results of azimuthal Fourier analyses of directed flow, vi vs. y/ncm [10,11]. 
Low energy (below 2 AGeV) data show no CM softening in neither or vi analyses. 
158 A.GeV Pb+Pb data show obvious and strong softening, ~ 100%. 
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